A simulation study of a windowless gas stripping room in an E//B neutral particle analyzer* 
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Neutral Particle Analyzer (NPA) is one of the crucial diagnostic devices on Tokamak facilities. Stripping unit 
is one of the main parts of the NPA. A windowless gas stripping room with two differential pipes is adopted in a 
parallel direction of electric and magnetic fields (E//B) NPA. The pressure distributions in the stripping chamber 
are simulated by Ansys Fluent together with MolFlow+. Based on the pressure distributions extracted from the 
simulation, the stripping efficiency of the E//B NPA is studied with GEANT4. The hadron reaction physics is 
modified to track the charge state of each particle in a cross section base method in GEANT4. The transmission 
rates (R) and the stripping efficiencies f+ı are examined for the particle energy ranging from 20 to 200 keV at 
the input pressure (Po) ranging from 20 to 400 Pa. According to the combined global efficiency, R x f+1, Po 
= 240 Pa is obtained as the optimum pressure for the maximum global efficiency in the incident energy range 


investigated. 
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I. INTRODUCTION 


Tokamak is a toroidal device used in nuclear fusion re- 
search for the magnetic confinement plasma. It provides 
a place to test the integrated technologies, materials, and 
physics regimes necessary for the future commercial produc- 
tion of fusion-based electricity [1]. Neutral Particle Analyzer 
(NPA) is one of the crucial diagnostic devices on Tokamak 
facilities. It is used to determine the bulk ion temperature, 
the isotopic ratio and the fast ion distribution of the plasma, 
by measuring the charge exchange neutral particles escaping 
from the plasma. Different types of NPA have been built 
in Tokamak facilities worldwide [2-13], such as the parallel 
direction of electric and magnetic fields (E//B) NPA on the 
Tokamak Fusion Test Reactor (TFTR) [2], the compact neu- 
tral particle analyzer (CNPA) on the Wendelstein 7-AS stel- 
larator [3], the low- and high- energy neutral particle analyz- 
ers (LENPA and HENPA) on the International Thermonuclear 
Experimental Reactor (ITER) [10], the solid state NPA (ss- 
NPA) on the Experimental Advanced Superconducting Toka- 
mak (EAST) [11] and the CP-NPA on the HuanLiugi-2A 
(HL-2A) [12]. 

Stripping unit plays an important role in the analyzing of 
the neutral particles, except for the flux measurement NPA 
such as ssNPA [11]. It provides a place to reionize the charge 
exchange neutral particles. According to the state of the strip- 
ping material, the stripping unit can be separated into two 
type, the stripping foil and the gas chamber. When a stripping 
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foil is used in the NPA for the low energy neutrals, an addi- 
tional accelerating or focusing voltage is required for the sec- 
ondary ions [3, 9, 10]. A carbon foil with the thickness of 100 
À is commonly used as the stripping foil. On the contrary, a 
gas chamber requires a differential pumping system when the 
stripping gas is used. Typically the integrated target thickness 
of the order of 1016 atoms/cm? for the H» gas is used in the 
Joint European Torus (JET) NPA [13], and 10!? atoms/cm? 
for the He gas is used in the E//B NPA on TFTR [2]. 

The energetic particles, also known as fast, superthermal, 
hot and high-energy particles, are expected to play a critical 
role in plasma heating, current drive, momentum transport, 
energy transfer and plasma stability [14, 15]. Many experi- 
mental and theoretical studies have been contributed to this 
field [16—20], and other related fields [21—26] recently. Aim- 
ing to study the frontier physics of the energetic particles and 
to measure the fuel ratio, a new E//B NPA has been designed 
in the present experimental devices [27]. This E//B NPA is a 
tandem type NPA like CNPA built in Ioffe Physicotechnical 
Institute, Russia [3]. It will provide mass resolution (H and D 
resolution) for particles in the energy range of 20 to 200 keV. 
The magnetic field is designed to be created with a permanent 
magnet for smaller size and simpler maintenance. The upper 
limit energy of the E//B NPA is determined from the negative 
ion source neutral beam heating on Huanliuqi-2M (HL-2M) 
device. The lower limit is set to 20 keV because we are in- 
terested in the fast ions, instead of the background ions. For 
more details, we refer to our previous work in Ref. [27]. 

In this article, the gas stripping chamber of the new E//B 
NPA is designed and studied. A windowless gas stripping 
chamber is adopted to avoid the replacement of the strip- 
ping foils and make an easy maintenance in the actual op- 
eration. The performance of the gas stripping chamber is in- 
vestigated using Ansys Fluent [28, 29] and MolFlow+ [30], 
together with GEANT4 [32, 33]. This article is organized as 
follows: The design and pressure calculation of the gas strip- 


ping chamber is presented in Sec. II. The results of GEANT4 
simulation and discussions are given in Sec. III. A brief sum- 
mary is given in Sec. IV. 


Il. DESIGN AND PRESSURE CALCULATION OF THE 
GAS STRIPPING CHAMBER 


Stripping unit is one of the main parts of the NPA. The 
electron of the escaped neutral particles will be stripped in 
the stripping unit. A windowless gas stripping chamber is 
adopted in the design to avoid the replacement of the strip- 
ping foil and for simple maintenance. In order to get a cer- 
tain high pressure inside the stripping room but high vacuum 
in the outside vacuum chamber at the same time, two differ- 
ential pipes with small flow conductance are used for strip- 
ping room. Fig. 1 shows the schematic layout of the stripping 
chamber. The stripping room (1) with two differential pipes 
(2) of 36 mm in length and 4 mm in diameter is placed inside 
a vacuum chamber (3). Two holes with a diameter of 6 mm (7 
and 8) are made on the entrance and exit flanges to limit the 
beam size and to maximize the vacuum isolation from the up- 
stream pipe and downstream chamber. H» gas is used as the 
stripping gas to avoid polluting the Tokamak fuel. It is filled 
in the stripping room from the top flange (5) of the vacuum 
chamber through a bellow (4). A machinery bearing molec- 
ular pump with a pumping speed of 340 L/s is used at the 
bottom of vacuum chamber together with a gate valve (11) in 
this work. 
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Fig. 1. Schematic diagram of the gas stripping chamber. 


The pressure distribution inside the stripping chamber is 
one of the main concerns of our design. The pressure 
of dozens of Pa is required in the stripping room to get 
enough stripping efficiency for high energy hydrogen (H) and 
deuteron (D) atoms. At this pressure region, the gas flow state 
in the stripping room stays at viscous-molecular flow [31]. 
The pressure in the gas inlet and the bellow is higher than 
that in the stripping room. On the contrary, two or three or- 
ders lower pressure is estimated in the vacuum chamber. The 


mean free path of the gas molecular inside the vacuum cham- 
ber is larger than the size of chamber and the motion of the 
gas molecular can be treated collisionlessly. The Monte Carlo 
code of MolFlow+ is often used to calculate the pressure dis- 
tribution of the collisionless gas in high vacuum system. But 
it is not accurate for all the gas region in the gas stripping 
chamber. The computational fluid dynamics (CFD) software, 
which includes the nonlinear effect of viscous fluid, has better 
performance in the high pressure region. However, the cal- 
culation of the CFD software at high vacuum region shows 
unphysical bump at corners. Therefore, the gas pressure dis- 
tribution of the stripping chamber is calculated, combining a 
CFD software of Ansys Fluent [28, 29] for the viscous re- 
gion in the bellow, stripping room, and differential pipes and 
a Monte Carlo software of MolFlow+ [30] for the low pres- 
sure collisionless region in the vacuum chamber. 


Three dimensional CFD calculations are performed in the 
Ansys Fluent software. The fluid region is established accord- 
ing to the structure of the stripping chamber shown in Fig. 1. 
The laminar vicious model is adopted in the calculations. A 
pressure type gas inlet is defined at the top flange (5) of vac- 
uum chamber and the pressure at the gas inlet Pp ranges from 
20 to 400 Pa with step of 20 Pa set in the simulation. Three 
pressure type gas outlets are defined at the entrance, exit holes 
(7 and 8) and the bottom of gate valve (11). A pressure of 
107 Pa is assumed for all three outlets. Due to the large flow 
conductance of the two differential pipes (2), small change in 
the outlet pressure does not affect the pressure distribution 
in the stripping room. Moreover, the pressure distribution 
in the vacuum chamber will be replaced with the results of 
MolFLow+. Therefore, the pressure of 10^? Pa of the outlets 
is used for all the Ansys Fluent calculations. The stainless 
steel is set as the wall material. The room temperature of 300 
K is used in the calculations. A typical gas flow rate of 9.97 
Pa-L/s is obtained at the bottom of the bellow for the input 
pressure Po = 100 Pa. 


In the low pressure region in the vacuum chamber (3) in 
Fig. 1 and inside the differential pipe (2) at |z| > 42 mm, 
where z = 0 is set to the central of the stripping room, the 
pressure distribution is simulated by MolFlow+ at the same 
temperature of 300 K. For the MolFlow- simulation, the out- 
gassing rate adopted is from the Ansys Fluent calculation at z 
= +42 mm inside the differential pipe, 4 mm from the pipe 
exit. For the pumping, it is assumed that the gas molecules 
are absorbed when they hit the surfaces of the entrance and 
exit holes (7 and 8), that is, the sticking factor is set to 1 on 
the surfaces. This results in the pumping speed of 12.4 L/s 
through the entrance and exit holes. A pumping speed of 340 
L/s is set at the bottom of the gate valve (11). 

The simulated two dimensional (2D) pressure distribution 
at Py = 100 Pa in the y-z plane at x = 0 is shown in Fig. 2 (a). 
A detailed pressure distribution around the stripping room is 
shown in Fig. 2 (b) in a magnified scale. The pressure distri- 
bution is presented along the beam line in Fig. 2 (c) plotted 
in a logarithmic scale. Using the two differential pipes in the 
design, the desired high pressure is achieved inside the strip- 
ping room, and the linear decreasing pressure is observed in- 
side the two differential pipes and the entrance and exit holes. 


A sharp change of the pressure at the entrance of the two dif- 
ferential pipes is observed in the Ansys Fluent calculations. 
It is not found when MolFlow- is used to simulate the whole 
gas region in the stripping chamber. 
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Fig. 2. (a) 2D pressure distribution in the y-z plane at x = 0 for the com- 
bined simulation of Ansys Fluent and MolFlow- at Po = 100 Pa. (b) Magni- 
fied pressure distribution around the stripping room. (c) Pressure distribution 
along the beam line as a function of z in a logarithmic scale. P1, P2, P3 and 
P4 indicate four typical pressures in the stripping chamber. 


To evaluate the pressure distribution changed as the pres- 
sure at the gas inlet changes, the pressures at four typical po- 
sitions are examined for all the Po investigated. Fig. 3 shows 
the pressure inside the stripping room (P), at entrance of the 
differential pipes (P2), at the vacuum chamber (3) and at 
the outside surfaces of entrance and exit holes (P4), which 
is pointed out in Fig. 2 (c), as a function of the pressure at 
the gas inlet P). Linear changes of Pj, P», P and Py on Po 
are obtained, though slight fluctuations are found. By using 
the two differential pipes in the design, more than 500 times 
lower pressure is observed in the vacuum chamber than that 
in the stripping room. Through the entrance and exit holes, 
about one order magnitude lower pressure is obtained for the 
upstream pipe and downstream chamber. 
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Fig. 3. Pressure P (solid circles), P» (solid squares), P3 (solid up triangles) 
and P4 (solid down triangles) as a function of the Po. The lines are obtained 
by the linear fits. 


The integrated target thickness (rr) is an important quan- 
tity of the gas stripping chamber. It is commonly used to 
evaluate the efficiency of a stripping chamber. Since the pres- 
sures at the four typical positions are used to construct the 
pressure distribution in GEANTA in the next section, a com- 


parison of the exact ny and the ny calculated from the pres- 
sure distribution used in GEANTA is necessary. Fig. 4 shows 
nr as a function of Po for the results of Ansys Fluent and 
MolFlow- (solid circles) and that of GEANT4 (open circles). 
Good agreement is found between them. 
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Fig. 4. n as a function of Po. Solid circles are the exact ny calculated 
from the pressure distributions of the Ansys Fluent and the MolFlow+ sim- 
ulations. Open circles are the ny calculated from the pressure distributions 
constructed in the GEANTA simulation. 


IIl. RESULTS OF GEANT4 SIMULATION AND 
DISCUSSIONS 


Utilizing the obtained Pj, P5, P3 and P^, the stripping 
efficiencies of the stripping room are further studied using 
GEANT4 [32, 33]. GEANT4 is a toolkit for the simulation 
of the passage of particles through matter. It has been applied 
for varies studies, such as our previous studies of the aver- 
age neutron detection efficiency for the DEtecteur MOdulaire 
de Neutrons detectors (DEMON) [35] and the module test 
of the Collision Centrality Detector Array (CCDA) [36], and 
for many other studies of the electron backscattering [34], the 
neutron time of flight spectrometer system at HL-2M [37], 
the performance of a large size CsI detector [38] and so on. 

In the GEANT4 simulation, the physics list includes the 
electromagnetic physics [39] and the hadronic physics [40], 
in which ion transportation, electromagnetic, nuclear elastic 
and inelastic processes are activated, though some of the pro- 
cesses may not be used in the actual simulation. Since the in- 
tegrated target thickness is around 1016 atoms/cm?, in which 
the scattering probability of the incident of particles and the 
target atoms is small enough that the multi-scattering is neg- 
ligible. Therefore, the G4ScreenedNuclearRecoil class [41, 
42] is included in the standard electromagnetic physics for 
the incident energy ranging from 10 eV to 100 MeV. 

The charge state of H and D atoms is the key variable in 
this study. However, the original GEANTA cannot handle the 
charge state properly. In order to simulate the charge state 
variation, the hadron reaction physics is modified to trace the 
charge state of each particle in a cross section base method. 
By introducing a global charge state variable in the hadron re- 
action physics, the charge state of H and D is recorded when 
the charge exchange reaction happened in the gas stripping 


chamber. Many charge exchange cross section measurements 
have been performed for H on Hə gas during the last cen- 
tury [43—53]. Fig. 5 shows the electron loss cross sections of 
H? (00,1), the electron capture cross sections of H ? (01,0), 
the electron capture cross sections of H® (00,—1) and the elec- 
tron loss cross sections of H~ (c..1,0) on He gas as a func- 
tion of the incident energy (E) in (a), (b), (c) and (d), respec- 
tively. Solid circles, solid squares, solid up triangles, solid 
down triangles, open circles, open squares and open up trian- 
gles represent the data of Gealy [43], Stier [44], Barnett [45], 
Sanders [46], Smith [47], McClure [48] and Van Zyl [49], re- 
spectively. Solid curves represent the ORNL recommended 
cross section [52, 53]. Since the experimental c... 9 does not 
cover the high energy region at E greater than 30 keV, the 
ORNL recommended cross section of H atom on H» gas is 
used in the simulation. For a given velocity, the charge ex- 
change cross sections of D on Cs [54] or Rb [55] vapor are 
the same as those of H. Approximately, the incident energy 
per nucleon (£/A) of H and D are the same for the energy 
range investigated when they have the same velocity. There- 
fore, the charge exchange cross sections of H are also used 
for D at the same E/A in the simulations. 
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Fig. 5. Charge exchange cross section of 09,1 (a), 71,0 (b), o0,-1 (c) 
and o_1,9 (d) in unit of 10-7 16cm? per molecule as a function of the in- 
cident energy (E) of H atom or ion. Solid circles, solid squares, solid up 
triangles, solid down triangles, open circles, open squares and open up trian- 
gles represent the data of Gealy [43], Stier [44], Barnett [45], Sanders [46], 
Smith [47], McClure [48] and Van Zyl [49], respectively. Solid curves rep- 
resent the ORNL recommended cross section [52, 53], which is also used in 
the simulation. 


The simulations are performed for H and D on H» gas with 
the incident energy ranging from 20 to 200 keV in the step 
of 20 keV, and with Po ranging from 20 to 400 Pa in the step 
of 20 Pa. H and D atoms are generated at the entrance hole 
(7) of the vacuum chamber, corresponding to the z position 
of -120 mm, and distributed uniformly on the entrance hole 
surface with the same diameter of 6 mm. The momentum 
direction is assumed parallel to the z-axis. One million events 
are generated for each run. The energy loss (AE/ E) of H and 
D at 20 keV as a function of P, is shown in Fig. 6. A slight 
lower energy loss is observed for D, because the mass of D is 


twice of that of H, which causes less energy loss during the 
collisions. One can see from Fig. 6 that the linear increasing 
trends are found for H and D as Po increases. The maximum 
energy loss of 20 keV H and D is less than 4% for all P; 
investigated. Comparing to the energy resolution of this NPA, 
the energy loss of H and D after passing through the stripping 
chamber is small and can be neglected. 
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Fig. 6. Energy loss of 20 keV H and D passing though the stripping chamber 
as a function of Po. Dashed line shows as the guide line of AE/E = 0.05. 
Solid lines are from the linear fits. 


Since the diameter of the two differential pipes is smaller 
than that of the entrance and the exit holes, some of the inci- 
dent particles will be stopped by the geometry of the stripping 
chamber. Moreover, due to the Coulomb scattering of the 
incident particles and the target atoms, some of the incident 
particles will be scattered away from their original directions. 
Therefore, the transmission rate (R) of the incident particles 
is important, especially for low energy particles which suffer 
more Coulomb scattering, when they pass though the strip- 
ping chamber. In this study, the transmission rate is defined 
as the ratio between the number of particles reached at the 
exit hole (8) after passing through the stripping chamber with 
the Hə gas (Po > 0 Pa) and that without the gas (Pp = 0 Pa, 
vacuum). In this way, the particle loss caused by the geom- 
etry of the stripping chamber is canceled out. Fig. 7 shows 
the transmission rate as a function of the incident energy for 
Po = 20 Pa (circles), 100 Pa (squares) and 400 Pa (triangles) 
in (a) and as a function of Po for the incident energy E = 20 
keV (circles), 40 keV (squares) and 100 keV (triangles) in (b). 
Solid and open symbols represent H and D, respectively. A 
slight increasing trend is observed for the transmission rate as 
the incident energy increases, but show an opposite trend as 
P increases. The scattering loss is small (less than 3%) for 
all the incident energies and the input pressures investigated. 

The stripping efficiency in the stripping chamber is evalu- 
ated using a charge fraction variable (f) for the incident par- 
ticles after the stripping area. The evolution of the charge 
fraction inside the stripping area is the most concerned in this 
study. Fig. 8 shows the charge state fraction as a function of 
Z position in the stripping chamber for H and D atoms at 20, 
100 and 200 keV. The pressures used in the simulation are 
those from the results of Po = 100 Pa as shown in Sec. II. 
Solid, dashed and long-dashed curves correspond to the frac- 
tion of the charge state 0, +1 and -1, respectively. Due to 
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Fig. 7. (a) R as a function of E. Circles, squares and triangles represent Po 
= 20, 100 and 400 Pa, respectively. (b) R as a function of Po. Circles, squares 
and triangles represent those at Po = 20, 40 and 100 keV, respectively. Solid 
and open symbols represent for H and D, respectively. Dashed lines show for 
FR = 1. The errors shown for the data points are the statistical errors. 


the small cross section of oo,.., the fraction of charge state 
-] is less than 2% for H and D at the incident energy of 20 
keV and becomes negligible for higher incident energies. The 
fraction of charge state 0 and +1 show a sharp change started 
from z around -40 mm, which is because the stripping room 
locates at -46 mm « z « 46 mm. Due to the energy depen- 
dence of the charge exchanging cross section, particles with 
lower incident energies have smaller saturation thickness for 
the charge fractions. A higher stripping efficiency for parti- 
cles with larger energies is mainly caused by the sharp de- 
crease of c o as the incident energy increases. 
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Fig. 8. Charge state fraction as a function of z position in the stripping 
chamber for H at 20 keV in (a), D at 20 keV in (b), H at 100 keV in (c), D at 
100 keV in (d), H at 200 keV in (e) and D at 200 keV in (f). The pressures 
used in the simulation are those from the results of Po = 100 Pa. Solid, 
dashed and long-dashed curves corresponding to fraction of charge state 0, 
+1 and -1, respectively. 


Figure 9 shows the fraction of charge state +1 (f41) atz = 
120 mm as a function of E/A for P) = 20 Pa (circles), 100 
Pa (squares) and 400 Pa (triangles) in (a), and as a function 
of Po for the incident energy E = 20 keV (circles), 100 keV 
(squares) and 200 keV (triangles) in (b). Solid and open sym- 
bols represent H and D, respectively. One can see from Fig. 9 
(a) that f41 increases at lower incident energies for different 


Po up to E/A around 100 keV. After reaching the maximum 
value, f,4 decreases for Py = 20 Pa, and stay flat for Po = 
100 Pa, but keep slowly increasing for Py = 400 Pa as E/A 
increases, which indicates that the thickness of the stripping 
gas is not enough for higher energy particles at Pp = 20 Pa. 
No noticeable difference between H and D is observed, indi- 
cating that the results can be applied also to neutral Tritium 
particles when E//A is used. As shown in Fig. 9 (b), the frac- 
tions of charge state +1 quickly reach a maximum value and 
keep the maximum as Pp increases at lower E. For larger EF, 
the fractions increase faster at lower Po and reach the maxi- 
mum at pressure around Po = 240 Pa. 
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Fig. 9. (a) f41 as a function of E/A. Circles, squares and triangles repre- 
sent the values at Po = 20, 100 and 400 Pa, respectively. (b) f+1 as a function 
of Po. Circles, squares and triangles represent those at Æ = 20, 100 and 200 
keV, respectively. Solid and open symbols represent H and D, respectively. 


In order to verify the GEANTA results, f; is also calcu- 
lated, using the gas integrated target thickness (nr) as 


fa- 


001 {1 


UMOR exp[-nmr(o01 + 010) - (1) 


The 691 and o1o are the stripping (the electron loss) and the 
charge exchange (the electron capture) cross sections of H(D) 
and H * (D*), respectively. The small amount of particle loss 
by the electron capture of H(D), o9_1 and o 10, is neglected 
in Eq. (1). The results are shown in Fig. 10. The symbols are 
the same f,, values from the GEANT4 simulation in Fig. 9 
(b) but plotted as a function of ny, and solid and dashed lines 
are those from Eq. (1) for H and D, respectively. Good agree- 
ments are found between the calculations and the GEANT4 
simulations. These good agreements originate from the fol- 
lowing fact. In the GEANT4 simulation, the evaluation of 
f41 is obtained from a Monte Carlo sampling of the charge 
state along the particle track according to the cross sections 
Of 091, 010. Co—1 and oio. As mentioned earlier, in Eq. (1) 
only a part of the cross sections (091 and c0) are used, ne- 
glecting o9_; and c. 10, since the latter values are orders of 
magnitude smaller. 

In order to determine the optimum condition, the global ef- 
ficiency, a combination of transmission rate and fraction of 
charge state +1, R x f4+1, is further studied. Fig. 11 shows 
the global efficiency of R x f,4 as a function of Po for the 
incident energy E = 20, 100 and 200 keV in (a), (b) and (c), 
and as a function of nr in (d), (e) and (f), respectively. The 
global efficiency of H and D decreases gradually as the pres- 
sure P) increases for E = 20 keV. On the other hand, the 
global efficiency shows similar trend as that of f,, as the 
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Fig. 10. The fraction of H+ (solid symbols) and D+ (open symbols) as 
a function of ny for E = 20 keV (circles), 100 keV (squares) and 200 keV 
(triangles). Solid and dashed lines are the results of Eq. (1). See more details 
in the text. 


pressure Po increases for E = 100 and 200 keV. For the pres- 
sure Py > 240 Pa, the global efficiency already becomes flat 
for all conditions of H and D at the incident energy E > 100 
keV. Considering the low temperature of the plasma in HL- 
2A/M, the number of high energy particles are in orders of 
magnitude less than that of the low energy particles. There- 
fore, Py = 240 Pa is obtained as the optimum pressure for the 
maximum global efficiency in the incident energy range in- 
vestigated. At this Po, the pressure at the vacuum chamber is 
less than 0.1 Pa, which is inside the operating pressure range 
of the molecular pump. The simulation results would provide 
a useful guide for the actual applications. 
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Fig. 11. Global efficiency of R x f+1 as a function of the pressure Po for 
the incident energy E = 20 keV (a), 100 keV (b) and 200 keV (c). (d), (e) 
and (f) are same as those of (a), (b) and (c) but as a function of ny. Solid 
and open circles represent for H and D, respectively. 


IV. SUMMARY 


Neutral particle analyzer (NPA) is one of the crucial di- 
agnostic devices on Tokamak facilities. The stripping unit 


is one of the main parts of NPA. A windowless gas tripping 
room with two differential pipes is adopted to keep a certain 
pressure for the parallel direction of electric and magnetic 
fields (E//B) NPA. The gas pressure distribution of the strip- 
ing chamber is calculated, combining a computational fluid 
dynamics (CFD) software of Ansys Fluent and a Monte Carlo 
software of MolFlow- for the low pressure collisionless re- 
gion in the vacuum chamber. The pressure distribution along 
the beam direction is obtained for different input pressures. 
A certain high pressure is achieved inside the stripping room 
and the linear decreasing pressure is obtained inside the dif- 
ferential pipes and the entrance and exit holes. More than two 
orders magnitude smaller pressure is obtained in the vacuum 
chamber than that inside the stripping room. 


Base on the pressure distributions calculated by Ansys Flu- 
ent and MolFlow-, the stripping efficiency of the stripping 
chamber for H and D atoms at the incident energy ranging 
from 20 to 200 keV is studied using GEANT4. The energy 
loss of H and D after passing through the stripping chamber 
is small and can be neglected for all the incident energies and 
the input pressures investigated. The scattering loss of H and 
D atoms on Hə gas is studied through the transmission rate 
(R) of the incident atoms. A slight increasing trend is ob- 
served for R as the incident energy increases, but show an 
opposite trend as the input pressure (Po) increases. The scat- 
tering loss is small (less than 3%) for all the incident energies 
and the input pressures investigated. 


A charge state variable is introduced to track the charge 
state of particles in the GEANT4 simulation. Adopting the 
ORNL recommended charge exchange cross sections in a 
modified hadron reaction physics, the charge state of each 
particle is traced in the simulation. The behaviors of charge 
fractions along the beam direction (z-axis) in the Hə gas are 
investigated for E = 20, 100 and 200 keV H and D atoms. The 
stripping efficiency is obtained as the fraction of charge state 
+1 at the exit hole of the vacuum chamber (z = 120 mm). Af- 
ter reaching the maximum value, f+; decreases for Pp = 20 
Pa, and stay flat for P = 100 Pa, but keep slowly increasing 
for Py = 400 Pa as the incident energy per nucleon increases. 
f+1 quickly reach a maximum value and keep the maximum 
as Po increases at lower FE. For larger F, the fractions in- 
crease faster at lower Po and reach the maximum at the input 
pressure around Po = 240 Pa. 


According to the combined global efficiency, R x f. Po 
— 240 Pa is found as the optimum pressure for the maximum 
global efficiency in the incident energy range investigated. 
The simulation results would provide a useful guide for the 
actual applications. 
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